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Abstract 

The most significant results on hard scattering processes in Au+Au, p+p, and 
d+Au collisions at ^snn = 200 GeV obtained after 3 years of operation at the 
BNL Relativistic Heavy-Ion Collider (RHIC) are summarized. Hadron pro- 
duction at high transverse momentum (pr) in central Au+Au collisions shows 
very different properties (production yields, flavor composition, azimuthal cor- 
relations, ...) in comparison to p+p, d+Au, and peripheral Au+Au collisions, 
and this fact constitutes arguably the most interesting outcome of the exper- 
imental program so far. Four main topics are covered in this report: (i) the 
depletion of high px inclusive hadron spectra, (ii) the "anomalous" baryon- 
meson composition at intermediate pr's, (iii) the disappearance of away-side 
azimuthal (jet) correlations, and (iv) the large value of collective azimuthal 
(elliptic flow) correlations. Additionally, the observed high pr "Cronin en- 
hancement" in d+Au collisions will be discussed. A succinct comparison of 
experimental data to theoretical models is also provided. 

1. Introduction 

Heavy-ion collisions at RHIC collider energies aim at the study of QCD matter at extreme energy den- 
sities. The driving force of this research program is the study of the transition from hadronic matter to 
a deconfined, chirally symmetric plasma of quarks and gluons (QGP) predicted to occur above a critical 
energy density e as 0.7 ± 0.3 GeV/fm 3 by lattice QCD calculations [1 1. Additionally, the combination 
of high-energies and large nuclear systems in the initial-state provides propitious grounds for the study 
of the regime of non-linear many-body parton (mostly gluon) dynamics at small-x where saturation ef- 
fects are expected to play a dominant role \2\. The most interesting phenomena observed so far in 
Au+Au reactions at ^s NN = 130 and 200 GeV are in the high transverse momentum (pr > 2 GeV/c) 
sector where the production of hadrons in central collisions shows substantial differences compared to 
more elementary reactions either in the "vacuum" (p+p, e + e~) or in a cold nuclear matter environment 
(d+Au). 

Several arguments make of hard scattering processes an excellent experimental probe in heavy-ion 
collisions: 

• They occur in the first instants of the reaction process (r ~ l/pr ^5 0.1 fm/c) in parton-parton 
scatterings with large Q 2 and, thus, constitute direct probes of the partonic phase(s) of the reaction. 

• A direct comparison to the baseline "vacuum" (p+p) production yields is straightforward after 
scaling by the nuclear geometry (in the simplest approach, this is done scaling by the nuclear 
overlap function Tab or, equivalently, by the number of nucleon-nucleon (NN) binary inelastic 
collisions, N co u = Tab • &nn, given by a Glauber model). 

• The production yields of high px particles are theoretically calculable via perturbative (or classical- 
field) QCD methods. 

Thus, high p? particles constitute experimentally and theoretically well calibrated observables that 
are sensitive to the properties of the dense QCD medium existing in central heavy-ion collisions. 

• Contribution to the CERN Yellow Report on Hard Probes in Heavy Ion Collisions at the LHC 
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The main observations so far at RHIC are the following: 

• The high pt yields of inclusive charged hadrons and 7r° in central Au+Au at y / sj^~ = 130 |3j|4ll3 
and 200 GeV 151171181191 1m . are suppressed by as much as a factor 4-5 compared to p+p and 
peripheral Au+Au yields scaled by Tab (or N co u). 

• At intermediate pr's (pt ~ 2. - 4. GeV/c) in central Au+Au, at variance with mesons (7r° HO 
and kaons fl3lD no suppression is seen for baryons (p,p ATI fT2l and A, A [13|), yielding an 
"anomalous" baryon over meson ratio p/ir ~ 1 much larger than the "perturbative" p/ir ~ 0. 1 - 
0.3 ratio observed in p+p collisions ff4lfT5ll and in e + e~ jet fragmentation [ 16 1. 

• The near-side azimuthal correlations of high p? (leading) hadrons emitted in central and peripheral 
Au+Au reactions [17 18] are, on the one hand, clearly reminiscent of jet-like parton fragmentation 
as found in p+p collisions. On the other, away-side azimuthal correlations (from back-to-back jets) 
in central Au+Au collisions are found to be significantly suppressed 1171 . 

• At low pt the strength of the azimuthal anisotropy parameter v 2 is found to be large and consistent 
with hydrodynamical expectations for elliptic flow. Above pr ~ 2 GeV/c where the contribu- 
tion from collective behaviour is negligible, V2 has still a sizeable value with a flat (or slightly 
decreasing) behaviour as a function of pt fT3l[T9l 20 1. 

• High pt production in "cold nuclear matter" as probed in d+Au reactions j2Tl l22l l23l ITO I not only 
is not suppressed but it is enhanced compared to p+p collisions, in a way very much reminiscent 
of the "Cronin enhancement" observed in p+A collisions at lower center-of-mass energies [24|. 

All these results point to strong medium effects at work in central Au+Au collisions, and have 
triggered extensive theoretical discussions based on perturbative or "classical"-neld QCD. Most of the 
studies on the high pt suppression are based on the prediction [ 25 1 that a deconfined and dense medium 
would induce multiple gluon radiation off the scattered partons, effectively leading to a depletion of the 
high-pr hadronic fragmentation products ("jet quenching"), though alternative interpretations have been 
also put forward based on initial-state gluon saturation effects ("Color Glass Condensate", CGC) 121 1261 . 
or final-state hadronic reinteractions [27]. The different behaviour of baryons and mesons at moderately 
high pt's has been interpreted, among others, in terms of "quark recombination" (or coalescence) effects 
in a thermalized partonic (QGP-like) medium |28], whereas the disappearance of the back-to-back 
azimuthal correlations can be explained in both QGP energy loss and CGC monojet scenarios. Finally, 
the large value of V2 above 2 GeV/c has been addressed by jet energy loss [29], gluon saturation [30|, 
and quark recombination I3T1 models. 

This summary report presents the px, ^snn, centrality, particle-species, and rapidity dependence 
of the inclusive high pr particle production, plus the characteristics of the produced jets and collective 
elliptic flow signals extracted from the azimuthal correlations at large px, as measured by the four exper- 
iments at RHIC (BRAHMS, PHENIX, PHOBOS and STAR) in Au+Au, p+p, and d+Au collisions. The 
whole set of experimental data puts strong constraints on the different proposed physical explanations 
for the underlying QCD medium produced in heavy-ion collisions at RHIC and at LHC energies. 



2. High pt hadron production in p+p collisions at \fs = 200 GeV 
2.1 p+p inclusive cross-sections: 

Proton-proton collisions are the baseline "vacuum" reference to which one compares the Au+Au results 
in order to extract information about the QCD medium properties. At y/s = 200 GeV, there currently exist 
three published measurements of high pt hadron cross-sections in p + p{p) collisions: UA1 p + p — > h^- 
(\rj\ < 2.5, pt < 7 GeV/c) [32 1, PHENIX p + p -> vr° (|r/| < 0.35, p T < 14 GeV/c) El, and STAR 
p + p -> h*- i\rj\ < 0.5, pt < 10 GeV/c) [7|. At yfs = 130 GeV, an interpolation between the ISR 
inclusive charged hadron cross-section and UA1 and FERMILAB data, has been also used as a reference 



for Au+Au at this value of y/s. Globally the spectra can be reasonably well parametrized by a power-law 
form A ■ (1 + px/po)~ n with the parameters 1 reported in Table l2.ll 
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Table 1: Parameters of the fit Ed' A a /dp 3 — A ■ (1 + pr/po) "to the inclusive pr distributions of all existing p + p(p) hadron 
(inelastic or non-singly diffractive) cross-sections measurements at y/s = 200 GeV. 



In general, all experimental results are consistent within each other, although it is claimed [ 7 1 that 
STAR p+p inclusive charged yield is smaller by a factor of 0.79 ±0.18 compared to UA1 p + p results 
(approximately independent of pr), the difference due in large part to differing non-singly-diffractive 2 
(NSD) cross section measured (35 ± 1 mb l3*2l in the first and 30.0 ± 3.5 mb [7 1 in the later). Standard 
next-to-leading-order (NLO) perturbative QCD calculations describe well the available high px p+p data 
at -/s = 200 GeV (see Fig. fflfor vr°). 
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Fig. 1 : High pr vr° cross-section in p+p collisions at y/s = 200 GeV (PHENIX) compared to the results of two different NLO 
pQCD calculations: (33) (left), (34) (right). 

2.2 p+p azimuthal correlations: 

PHENIX [ 35 | has studied the azimuthal correlations at high px in p+p collisions at y/s = 200 GeV 
extracting several parameters characterizing the produced jets: 

• Mean jet fragmentation transverse momentum: = 373 ± 16 MeV/c, in agreement with 

previous measurements at ISR [36] and showing no significant trend with increasing yfs. 

'The fit parameters po and n are actually strongly correlated via the mean pr of the collision: (pr) = 2po /(n — 3). 
2 PHENIX high pr 7r° cross-section is inclusive and contains, in principle, all inelastic (including diffractive) channels. 



• Average partem transverse momentum (fitted to a constant above 1.5 GeV/c): (\k± y \) = 725 
± 34 MeV/c. The momentum of the pair p± is related to the individual parton (IfeijJ) via 

J {\p\\) P air = (\k± y \). The extracted yj (\p]_\)pair = 1.82 ± 0.85 GeV/c is in agreement 
with the existing systematics of dimuon, diphoton and dijet data in hadronic collisions [37 ]. 



3. High pt hadron production yields in Au+Au collisions 

There is a significant amount of high px Au+Au experimental spectra (px > 2 GeV/c) measured by the 
4 experiments at RHIC: inclusive charged hadrons at 130 011113 and 200 GeV (7J HI El EDI, neutral 
pions at 130 and 200 GeV J6), protons and antiprotons at 130 ifTTI and 200 GeV [12), K° s at 200 
GeV 1131 . and A, A at 200 GeV fT3l . Moreover, all these spectra are measured for different centrality 
bins and permit to address the impact parameter dependence of high px production. 

Details on hadron production mechanisms in AA are usually studied via their scaling behavior 
with respect to p+p collisions. On the one hand, soft processes (px < 1 GeV/c) are expected to scale 
with the number of participating nucleons N part l38| (and they actually approximately do [39 40 1). On 
the other, in the framework of collinear factorization, hard processes are incoherent and thus expected to 
scale with N co n. 

The first interesting result at RHIC in the high px sector is the breakdown of this N co u scaling for 
central Au+Au collisions. Fig. [2] shows the comparison of the p+p ir° spectrum to peripheral (left) and 
central (right) Au+Au spectra, and to pQCD calculations. Whereas peripheral data is consistent with a 
simple superposition of individual AW collisions, central data shows a suppression factor of 4 - 5 with 
respect to this expectation. 




Fig. 2: Invariant n yields measured by PHENIX in peripheral (left) and in central (right) Au+Au collisions (stars), compared 
to the Ncou scaled p+p tv° yields (circles) and to a NLO pQCD calculation |41 1 (gray line). The yellow band around the scaled 
p+p points includes in quadrature the absolute normalization errors in the p+p and Au+Au spectra as well as the uncertainties 
in Tab ■ Updated version of Fig. 1 of 1 42 1 with final published data 1 511331. 



It is customary to quantify the medium effects at high px using the nuclear modification factor 
given by the ratio of the A A to the p+p invariant yields scaled by the nuclear geometry (Tab)' 

d?Nf A /dydp T 



{Tab) x d 2 af /dydp 



T 



Raa{pt) measures the deviation of AA from an incoherent superposition of NN collisions in terms of 
suppression (Raa <1) or enhancement (Raa >!)• 



3.1 High pt suppression: magnitude and px dependence 

Figure |3] shows Raa(pt) for (STAR Q, left) and tt° (PHENIX [6|, right) measured in peripheral 
(upper points) and central (lower points) Au+Au reactions at ^s^n = 200 GeV. As seen in Fig. |2j 
peripheral collisions are consistent 3 with p+p collisions plus N co u scaling as well as with standard pQCD 
calculations [ 43 , 44 1 , while central Au+Au are clearly suppressed by a factor ~ 4 - 5. 
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Fig. 3: Nuclear modification factor, Raa(pt), in peripheral and central Au+Au reactions for charged hadrons (left) and 7T° 
(right) measured at ^/sjvjv = 200 GeV by STAR and PHENIX respectively. A comparison to theoretical curves: pQCD-I|43|, 
pQCD-II |44], is also shown. 

The high px suppression in central collisions for both 7r° and h is smallest at px = 2 GeV/c and 
increases to an approximately constant suppression factor of I/Raa ~ 4 - 5 over px = 5 - 10 GeV/c. 
Above 5 GeV/c the data are consistent within in errors with "participant scaling" given by the dotted line 
at Raa ~ 0.17 in both plots (actually both STAR and PHENIX data are systematically slightly above 
this scaling). The magnitude and pt dependence of Raa in the range px = 1 - 10 GeV/c (corresponding 
to partem fractional momenta x\p = Pt I ' \fs{e. ±Vl + e ±m ) « 2px/ ' yfs ~ 0.02 - 0.1 at midrapidity), 
is alone inconsistent with "conventional" nuclear effects like leading-twist shadowing of the nuclear 
parton distribution functions (PDFs) 1481 l49l . Different pQCD-based jet quenching calculations B3l 
El 113 HU EID based on medium-induced radiative energy loss, can reproduce the magnitude of the tt° 
suppression assuming the formation of a hot and dense partonic system characterized by different, but 
related, properties: i) large initial gluon densities dN 9 /dy ~ 1000 [43 1, ii) large "transport coefficients" 
% ~ 3.5 GeV/fm 2 [45 |, hi) high opacities L/X ~ 3 - 4 ll34l . or iv) effective parton energy losses of the 
order of dE/dx ~ 14 GeV/fm ll44l . 

The px dependence of the quenching predicted by all models that include the QCD version of 
the Landau-Pomeranchuck-Migdal (LPM) interference effect (BDMPS [50 1 and GLV [51 1 approaches) 
is a slowly (logarithmic) increasing function of px, a trend not compatible with the data over the entire 
measured px range. Other approaches, such as constant energy loss per parton scattering, are also not 
supported as discussed in BT1 . Analyses which combine LPM jet quenching together with shadowing 
and initial-state px broadening ("pQCD-II" [43 1 in Fig. |3jl globally reproduce the observed flat px 
dependence of Raa, as do recent approaches that take into account detailed balance between parton 
emission and absorption ("pQCD-I" [44 1 in Fig. |3] left). 

At variance with parton energy loss descriptions, a gluon saturation calculation [26 1 is able to 
predict the magnitude of the observed suppression, although it fails to reproduce exactly the flat px 

3 Although peripheral STAR charged hadron data seems to be slightly above Raa = 1 and PHENIX ir° data seems to be 
below, within errors both measurements are consistent with "collision scaling". 



dependence of the quenching [7]. Similarly, semi-quantitative estimates of final-state interactions in a 
dense hadronic medium *2*71 yield the same amount of quenching as models based on partonic energy 
loss, however it is not yet clear whether the pt evolution of the hadronic quenching factor is consistent 
with the data or not f71l53ll. 

The amount of suppression for vr° [6] and |7 8] is the same above px ~ 4 - 5 GeV/c for all 
centrality classes 0. However, below px ~ 5 GeV/c, 7r°'s are more suppressed than inclusive charged 
hadrons in central collisions (as can be seen by comparing the right and left plots of Fig. |3j). This is due 
to the enhanced baryon production contributing to the total charged hadron yield in the intermediate pt 
region (pr ~ 1 - 4 GeV/c) in Au+Au collisions [12 13 ] (see section l3~4l below'). 



3.2 High pt suppression: ^/s NN dependence 

Figure |4] shows Raa{pt) for several ir° measurements in high-energy AA collisions at different center- 
of-mass energies f52l . The PHENIX Raa(pt) values for central Au+Au collisions at 200 GeV (circles) 
and 130 GeV (triangles) are noticeably below unity in contrast to the enhanced production (Raa >1) 
observed at CERN-ISR (min. bias a + a 031, stars) and CERN-SPS energies (central Pb+Pb 14*71 . 
squares) and interpreted in terms of initial-state p? broadening ("Cronin effect" [24 1). 



< 
< 



2.5 

2 
1.5 

Y 
0.5 



7t° 0-10% central: 






• Au+Au @ 




= 200 GeV 


Au+Au @ 


*J S NN 


= 130 GeV 


■fc a + a @ * 


J S NN 


= 31.0 GeV 


■ Pb+Pb @ 




= 17.3 GeV 




8 10 
p T (GeV/c) 



Fig. 4: Nuclear modification factor, Raa(pt), for n° measured in central ion-ion reactions at CERN-SPS 1471 (squares), 
CERN-ISR [15 ! (stars), and BNL-RHIC (triangles |3 1, circles Q>)) energies. 



Figure [5] shows R2oo/13o{pt), the ratio of Au+Au charged hadron yields at y / s^7 =130 and 
200 GeV in 4 centrality classes, compared to pQCD and gluon saturation model predictions [7|. The 
increase in high p-r yields between the two center-of-mass energies is a factor ~2 at the highest pr's, 
whereas at low pr, the increase is much moderate, of the order of 15%. The large increment of the hard 
cross sections is naturally consistent with pQCD expectations due to the increased jet contributions at 
high transverse momenta. In the saturation model [26 1 the increase at high pr is accounted for by the 
enhanced gluon densities at ^/s NN = 200 GeV compared to 130 GeV in the "anomalous dimension" xt 
region of the Au parton distribution function. 

PHENIX [ 8 ] has addressed the yfs dependence of high p? production by testing the validity of 
"xt scaling" in Au+Au, i.e. verifying the parton model prediction that hard scattering cross sections can 
be factorized in 2 terms depending on yfs and xt = %Pt/ 'yfs respectively: 

^XT = ^V F (^) ^ E ^ = ?-^G(x T ). (2) 

dp p™(v«) dp 3 ^n(x T ,Vs) 

In Q, F(xt) embodies all the xt dependence coming from the parton distribution (PDF) and frag- 
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Fig. 5: R20o/13o{pt) vs. pr for (h + + h )/2 for four different centrality bins measured by STAR compared to pQCD and 
gluon saturation model predictions |7 1. 



mentation (FF) functions 4 , while the exponent n, related to the underlying parton-parton scattering, is 
measured to be n « 4 - 8 in a wide range of p + p, p collisions [8]. Fig.[6]compares the xy-scaled hadron 
yields in ^snn = 130 GeV and 200 GeV Au+Au central and peripheral collisions. According to Eq. 
the ratio of inclusive cross-sections at fixed xy should equal (200/130)". On the one hand, xt scaling 
holds 5 in Au+Au with the same scaling power n = 6.3 ± 0.6 for neutral pions (in central and peripheral 
collisions) and charged hadrons (in peripheral collisions) as measured in p+p [8 1. This is consistent with 
equal (pQCD-like) production dynamics in p+p and Au+Au, and disfavours final-state effects described 
with medium-modified FF's that violate xt scaling (e.g. constant partem energy losses independent of 
the parton pr). Equivalently, models that predict strong initial-state effects (e.g. gluon saturation) respect 
xt scaling as long as their predicted modified nuclear PDFs are depleted, independently of \/s, by the 
same amount at a given xt (and centrality). On the other hand, Fig. |6] (right) shows that charged hadrons 
in central collisions (triangles) break xt scaling which is indicative of a non perturbative modification of 
particle composition spectra from that of p+p at intermediate pr's (see section l3~4l below). 
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and 200 GeV by PHENIX 1 8 1. Central (Peripheral) xt spectra are represented by triangles (squares), and solid (open) symbols 
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130 GeV scaled by a factor of [130/200] 6 3 ). 



4 PDFs and FFs, to first order, scale as the ratio of pr at different yfs. 
5 In the kinematical region, xt > 0.03, where pQCD is expected to hold. 



3.3 High pt suppression: centrality dependence 

In each centrality bin, the value of the high px suppression can be quantified by the ratio of Au+Au 
over A^co^-scaled p+p yields integrated above a given (large enough) px- The centrality dependence 
of the high px suppression for 7r° and charged hadrons, given by Raa(pt > 4.5 GeV/c), is shown in 
Fig. (left) as a function of (N paTt ) for PHENIX data. The transition from the N co u scaling behaviour 
(Raa ~ 1) apparent in the most peripheral region, (N part } < 40, to the strong suppression seen in central 
reactions (Raa ~ 0.2) is smooth. Whether there is an abrupt or gradual departure from N co u scaling in 
the peripheral range cannot be ascertained within the present experimental uncertainties [52]. The data, 
however, is inconsistent with N co u scaling (at a 2a level) for the 40-60% centrality corresponding to 
(N par t} ~ 40 - 80 (5l , whose estimated "Bjorken" energy density (esj ~ 1 GeV/fm 3 ) 1521 is in the 
ball-park of the expected "critical" QCD energy density. A similar centrality dependence of the high px 
suppression is seen m STAR h± data (Fig.El right) 
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Fig. 7: Evolution of the high pr ir° and h suppression, Raa{pt > 4.5 GeV/c), as a function of centrality given by {N par t) 
(PHENIX, left). Same evolution shown as Raa{pt >6.0 GeV/c) for STAR h ± data (53) (right). 



Npart (instead of N co u) scaling at high px is expected in scenarios dominated either by gluon 
saturation [26| or by surface emission of the quenched jets ll54l . "Approximate" N par t scaling has been 
claimed by PHOBOS [ 9 1 : the ratio of central to a fit to mid-central yields in the range px ~ 2. - 4. GeV/c 
stays flat as a function of centrality (Fig. [S] left). However, at higher px values, where the suppression 
is seen to saturate at its maximum value, the centrality dependence of the ratio of A^a^-scaled Au+Au 
over p+p yields for 7r° and measured by PHENIX (Fig. |8j right) does not show a true participant 
scaling (R P aa > 1 f° r an centralities). Nonetheless, the fact that the production per participant pair 
above 4.5 GeV/c is, within errors, approximately constant over a wide range of intermediate centralities, 
is in qualitative agreement with a gluon saturation model prediction \26\. 



3.4 High px suppression: particle species dependence 

One of the most intriguing results of the RHIC program so far is the different suppression pattern of 
baryons and mesons at moderately high pr's. Figure |9] (left) compares the N co u scaled central to pe- 
ripheral yield ratios 6 for (p + p)/2 and vr°: R cp = (yiel<p- w% ^ /iV c ° o 7/ 0% ) /(yielS 60 - 92 ^ /N™J 92% ). 
From 1.5 to 4.5 GeV/c the (anti)protons are not suppressed (R cp ~ 1) at variance with the pions which 
are reduced by a factor of 2 - 3 in this px range. If both 7r° and p, p originate from the fragmentation 
of hard-scattered partons that lose energy in the medium, the nuclear modification factor R cp should 
be independent of particle species contrary to the experimental result. The same discussion applies for 
strange mesons and baryons as can be seen from the right plot of Fig. [9j Whereas the kaon yields in 

6 Since the 60-92% peripheral Au+Au (inclusive and identified) spectra scale with iVcoii when compared to the p+p yields 1 6 
I7I I12I . R cp carries basically the same information as Raa- 
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Fig. 8: Left: Ratio of central to semi-central h yields normalized by N par t in 4 different pr bins, as a function of N part 
measured by PHOBOS [9 1. The dashed (solid) line shows the expectation for N co u (Npart) scaling. Right: Ratio of Au+Au 
over p+p 7T° and h ± yields above 4.5 GeV/c normalized by N par t (Raa™) as a function of centrality given by {N par t) as 
measured by PHENIX 1 6 8 1. The dashed line indicates the expectation for N par t scaling. 



central collisions are suppressed with respect to "N co u scaling" for all measured pt, the yield of A + A 
is close to expectations from collision scaling in the p? range 1.8-3.5 GeV/c. Interestingly, above pj> ~ 
5.0 GeV/c, the K s , K^, A + A, and charged hadron yields are suppressed from binary scaling by a 
similar factor. 
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Au+Au collisions: (p + p)/2 (dots) and tv° (circles) by PHENIX [HI (left), and A, A (circles) and K°,K ± (triangles) by 
STAR EH (right). 



Figure^] (left) shows the ratios of (p + p)/2 over 7r° as a function of px measured by PHENIX 
in central (0-10%, circles), mid-central (20-30%, squares), and peripheral (60-92%, triangles) Au+Au 
collisions [12], together with the corresponding ratios measured in p+p collisions at CERN-ISR ener- 
gies I I4IIT5II (crosses) and in gluon and quark jet fragmentation from e + e~ collisions [ 16 1 (dashed and 
solid lines resp.). Within errors, peripheral Au+Au results are compatible with the p+p and e + e~ ra- 
tios, but central Au+Au collisions have a p/ir ratio ~ 4 - 5 times larger. Such a result is at odds with 
standard perturbative production mechanisms, since in this case the particle ratios p/it and p/n should 



be described by a universal fragmentation function independent of the colliding system, which favors 
the production of the lightest particle. Beyond px ~ 4.5 GeV/c, the identification of charged particles 
is not yet possible with the current PHENIX configuration, however the measured h/ir° ~ 1.6 ratio 
above px ~ 5 GeV/c in central and peripheral Au+Au is consistent with that measured in p+p collisions 
(Fig. QUI right). This result together with STAR R cp result on strange hadrons (Fig. |9j right) supports 
the fact that for large px values the properties of the baryon production mechanisms approach the (sup- 
pressed) meson scaling, thus limiting the observed baryon enhancement in central Au+Au collisions to 
the intermediate transverse momenta px < 5 GeV/c. 

Several theoretical explanations (see refs. in fi"Tl[T2"ln have been proposed to justify the differ- 
ent behaviour of mesons and baryons at intermediate p^'s based on: (i) quark recombination (or coa- 
lescence), (ii) medium-induced difference in the formation time of baryons and mesons, (iii) different 
"Cronin enhancement" for protons and pions, or (iv) "baryon junctions". In the recombination picture 
[ 28 1 the partons from a thermalized system coalesce and with the addition of quark momenta, the soft 
production of baryons extends to much larger values of px than that for mesons. In this scenario, the 
effect is limited to px < 5 GeV, beyond which fragmentation becomes the dominant production mecha- 
nism for all species. 
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3.5 High px suppression: pseudorapidity dependence 

BRAHMS is, so far, the only experiment at RHIC that has measured high px inclusive charged hadron 
spectra off mid-rapidity. Fig. ^J(left) shows the nuclear modification factors Raa(pt) for central and 
semi -peripheral Au+Au measurements at mid-pseudorapidity (rj = 0) and at rj = 2.2 [ 10 1 . The high px 
suppression is not limited to central rapidities but it is clearly apparent at forward 77's too. Fig. ITTl(right) 
shows the ratio of suppressions at the two pseudorapidity values, R v = R cp (rj = 2.2) /R cp (r) = 0). The 
high px deficit at 77 = 2.2 is similar to, or even larger than, at r\ = indicating that the volume causing the 
suppression extends also in the longitudinal direction. 

4. High px azimuthal correlations in Au+Au collisions 

There are two main sources of azimuthal correlations at high px in heavy-ion collisions: 

• The fragmentation of hard-scattered partons results in jets of high px hadrons correlated in both 
rapidity and azimuthal angle. Such correlations are short range (A77 < 0.7, < 0.75), involve 




Fig. 1 1 : Left: Raa (pt) measured by BRAHMS at 77 = and 77 = 2.2 for 0-10% most central and for semi-peripheral (40-60%) 
Au+Au collisions. Right: Ratio R v of R cp distributions at 77 = 2.2 and 77 = 0. Figs, from 1101 . 



comparatively large transverse momentum particles (pr > 2 GeV/c), and are unrelated (in princi- 
ple) to the orientation of the AA reaction plane. 
• Collective (elliptic) flow: The combination of (i) the geometrical asymmetry in non-central AA re- 
actions ("almond"-like region of the overlapping nuclei), and (ii) multiple reinteractions between 
the produced particles in the overlap region; results in pressure gradients in the collision ellipsoid 
which transform the original coordinate-space asymmetry into a momentum-space anisotropy. The 
amount of elliptic flow (a true collective effect absent in p+p collisions) is measured by the sec- 
ond harmonic coefficient, V2 = (cos(20)), of the Fourier expansion of the particles azimuthal 
distribution with respect to the reaction plane. 

Additionally, there are other second-order sources of angular correlations like resonance decays, 
final state (particularly Coulomb) interactions, momentum conservation, or other experimental effects 
like photon conversions, which have to be subtracted out in order to extract the interesting "jet-like" or 
"flow-like" signals. 



4.1 High pt azimuthal correlations: Jet signals 

Although, standard jet reconstruction algorithms fail below pt ~ 40 GeV/c when applied to the soft- 
background dominated environment of heavy-ion collisions, angular correlations of pairs of high pt 
particles [17 18 1 have been very successfully used to study on a statistical basis the properties of the 
produced jets. For each event with "trigger" particle(s) with p T = 4 - 6 GeV/c and "associated" particle(s) 
with px = 2 - 4 GeV/c and |r/| < 0.7, STAR [ 17 1 determines the two-particle azimuthal distribution 

D(Acp) oc - 1 f" . (3) 

^trigger d{A(p) 

Fig. [\2~\ shows D(A(j)) for peripheral (left) and central (right) Au+Au collisions (dots) compared to 
D(A(p) from p+p collisions (histogram), and to a superposed cos(2A<ft) flow-like term (blue curve). 
On the one hand, the correlation strength at small relative angles (A<p ~ 0) in peripheral and central 
Au+Au as well as at back-to-back angles (A^ ~ it) in peripheral Au+Au are very similar to the scaled 
correlations in p+p collisions. The near-side peaks in all three collision systems are characteristic of jet 
fragmentation [ 17 ] (a result also observed by PHENIX using neutral trigger particles [ 18 1). On the other 
hand, the away-side peak (Acf) ~ tt) in central collisions is strongly suppressed. 

In order to study the evolution as a function of centrality of the the near-side, D u (A(p < 0.75), 
and away-side, D AuAu (A<ft > 2.24), angular correlations in Au+Au compared to p+p, D pp , STAR has 




Fig. 12: Azimuthal correlations for peripheral (left) and central (right) Au+Au collisions compared to the pedestal and flow- 
scaled correlations in p+p collisions. Fig. from 1551 . 



constructed the quantity 



jgg d{^)[D AuAu ~ BQ. + 2j co S (2A^))] 



(4) 



where B accounts for overall background and V2 the azimuthal correlations due to elliptic flow. Fig. ^] 
shows Iaa for the near-side (squares) and away-side (circles) correlations as a function of the number 
of participating nucleons (N part ). On the one hand, the near-side correlation function is relatively sup- 
pressed compared to the expectation from Eq. © in the most peripheral region (a result not completely 
understood so far) and increases slowly with N paTt . On the other hand, the back-to-back correlation 
strength above the background from elliptic flow, decreases with increasing N part and is consistent with 
zero for the most central collisions. The disappearance of back-to-back jet-like correlations is consistent 
with large energy loss effects in a system that is opaque to the propagation of high-momentum partons 
or their fragmentation products. 
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4.2 High pt azimuthal correlations: Collective elliptic flow 

At low pt the strength of the elliptic flow signal is found to be large and consistent with hydrodynam- 
ics expectations. Above pr ~ 2 GeV/c where the contribution from collective behaviour is negligi- 
ble, V2 is found to be still a sizeable signal which saturates and/or slightly decreases as a function of 
Pt fT9l l20l IT3l . The large value V2{pt > 2 GeV/c) ~ 0.15 implies unrealistically large partem densities 
and/or cross-sections according to standard parton transport calculations ll3Tl . Various interpretations 
have been proposed to account for such a large V2 parameter within different physical scenarios. In jet 
quenching models [29 ] the resulting momentum anisotropy results from the almond-like density pro- 
file of the opaque medium (see, however, l56l ). Calculations based on gluon saturation 11301 yield a 
("non-flow") azimuthal asymmetry component from the fragmentation of the released gluons from the 
initial-state saturated wave functions of the colliding nuclei. Finally, quark recombination effects l3"Tl can 
naturally enhance the elliptic flow of the produced hadrons compared to that of partons. The measured 
V2(pt) for mesons and baryons shows a distinct pattern (Fig.[T4li: v™ > v\ at low pr, v™ ~ v\ at pt ~ 
2 GeV/c, and v™ < v\ at higher pt's; which further constraints the proposed theoretical explanations. 



0.3 



0.2 



0.1 



' I ' I ' I ' 

° K~ 

• pbar -^j 


-* P - 


r^Pi 






m 
E 
2 

ro 0.2 



B 01 
o 
w 
'c 

< 





° A + A 


o h + +h~ 




Hydro calculations 








K 








K 

~ - - - P 
A 






















.•'.a 

















2 3 4 
p T (GeV/c) 



2 3 4 
p T (GeV/c) 



2 4 

Transverse Momentum p T (GeV/c) 



Fig. 14: V2 as a function of transverse momentum for identified particles at RHIC: ir , K and p, p from PHENIX (left), and 
K° and A, A from STAR (right). 

Quark coalescence models [31 ] naturally lead to weaker baryon flow than meson flow at low pt, 
while the opposite holds above 2 GeV/c. This simple mass ordering expectation is confirmed by the 
identified particle data from PHENIX and STAR (Fig. [T31 . The fact that the V2 parameters scaled by the 
number of constituent quarks (n = 2 for mesons, n = 3 baryons) versus pT/n, globally fall in a single 
curve, supports the scenario where hadrons at moderate pt's form by coalescence of co-moving quarks. 
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5. High pt hadron production in d+Au collisions 

Proton- (or deuteron-) nucleus collisions constitute a reference "control" experiment needed to determine 
the influence of cold nuclear matter effects in high pt hadro-production. Since final-state medium effects 
are marginal in p,d+Au collisions, they are basic tools to ascertain whether models based on initial- or 
final- state QCD effects can explain the distinct hard scattering behaviour observed in Au+Au collisions 
at RHIC. During the third year of RHIC operation, the 4 experiments collected data from d+Au collisions 
at ^/s^ = 200 GeV. The resulting high p T results at mid-rapidity from PHENIX ED, STAR E3, 
PHOBOS [23|, and BRAHMS 033 consistently indicate the following: 

• High px inclusive I2T1 l22l l23l fTOl and 7r° [21 1 spectra from d+Au minimum bias (MB) col- 
lisions are not suppressed but are enhanced compared to p+p collisions (RdAu plots in Fig. fl6l> . 
in a way very much reminiscent of the "Cronin effect" observed in fixed-target p+A collisions at 
lower yfs [24]. As a matter of fact, p+Au collisions (from neutron-tagged d+Au events I2T1D show 
a similar behaviour as minimum bias d+Au collisions. 

• Above pt ~ 2.5 GeV/c the nuclear modification factor of inclusive charged hadrons in MB d+Au 
collisions saturates at 7 R^a 1.4. Above 6 GeV/c, STAR /i ± and PHENIX vr° results seem to 
indicate that RdAu decreases as a function of pt, becoming consistent with 1 at around 8 GeV/c. 

• The "Cronin enhancement" for unidentified hadrons at high pt (RdAu ra 1-35) is larger than for 
neutral pions (Rf Au « 1.1) ED- 

• The degree of "enhancement" in d+Au compared to p+p increases with collision centrality l23l 
22], an opposite trend to Au+Au results. 

• The azimuthal correlations in MB and central d+Au collisions are very similar to that of p+p 
and do not show the significant suppression of the away-side peak observed in central Au+Au 
reactions 1221 . 

All these results lead to the conclusion that no "cold" nuclear matter (or initial-state) effects, - 
like a strong saturation of the nuclear parton distribution functions in the relevant (x, Q 2 ) kinematical 
region probed by the current experimental setups-, can explain the high pt behaviour in central Au+Au. 
The data suggest, instead, that final-state interactions are responsible of the high pt suppression and the 
disappearance of back-to-back dijet correlations observed at mid-rapidity in central Au+Au reactions. 

6. Summary 

A vast body of high transverse momentum data has been collected and analyzed during 2000-2003 by the 
4 experiments at RHIC in Au+Au, p+p, d+Au collisions at Js NN = (130)200 GeV. This has permitted 
a detailed comparative study of the properties of high pt particle production in high-energy heavy-ion 
collisions as a function of pt, ^/s NN , collision centrality, pseudo-rapidity, and particle species. All these 
studies reveal that hadron production at high transverse momentum in central Au+Au reactions shows 
significant deviations compared to p+p, d+Au, and Au+Au peripheral collisions at RHIC energies, as 
well as to nucleus-nucleus data at lower center-of-mass energies. The main observations are: (i) the 
inclusive and identified spectra are suppressed by a factor 4-5 above pt ~ 5 GeV/c compared to the 
expectations of N co u scaling, (ii) the baryon yields are enhanced compared to the meson yields in the 
range of intermediate momenta pt ~ 2 - 5 GeV/c, (iii) the back-to-back dijet azimuthal correlations are 
significantly suppressed, (iv) there is a strong constant elliptic flow signal at high pr- The whole set of 
results puts strong experimental constraints on the properties of the underlying QCD medium produced 
in Au+Au reactions at collider energies. Comparison of the energy spectra and angular correlations data 
to the theoretical calculations globally supports pQCD-based models of final-state parton energy loss in 
a dense medium, although non-perturbative effects (like e.g. quark coalescence) are needed in order to 
explain the baryon-meson differences in yield and V2 in the intermediate pt window (j>t ~ 2 - 5 GeV/c). 

7 KTu 1X (Pt = 2-7 GeV/c) ~ 1.35, R^ipr = 2-6 GeV/c) ~ 1.45, #^ hms (pt = 2-5 GeV/c) ~ 1.3. 
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Fig. 16: Top: Nuclear modification factor R dAu for MB d+Au: and tt° (PHENIX, Ze/f), ft* (BRAHMS, n'gfo). Bot- 
tom: -RdAu for ft measured by PHOBOS in 4 different d+Au centralities (left), and comparison of two-particle azimuthal 
distributions for central d+Au, p+p, and central Au+Au collisions (STAR, right). 

Theoretical predictions of a strong saturation of the nuclear wave function at high energies are also in 
agreement with most of the data but do not seem to explain consistently Au+Au and d+Au RHIC results 
at midrapidity. Coming ion-ion runs at RHIC and, in the mid-term, Pb+Pb collisions at LHC energies 
will help to further strengthen our current understanding of the physics of QCD media at high energy 
densities. 
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